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Shape anisotropy provides a simple mechanism to adjust the local bias field in pat- 
terned structures. It is well known that for ellipsoidal particles < 1 fim in size, a 
quasi-single domain state can be realized with uniform anisotropy field. For larger 
patterned ferromagnetic thin-film elements, domain formation is thought to limit the 
effectiveness of shape anisotropy. In our work we show that very soft lithographically 
patterned Ni8oFe2o films with control of induced magnetic anisotropy can exhibit 
shape anisotropy fields in agreement with single-domain models, for both hysteresis 
loop measurements at low field and ferromagnetic resonance measurements at high 
field. We show the superiority of the fiuxmetric form over the magnetometric form 
of anisotropy estimate for thin films with control dimensions from 10 /im to 150 fim 
and in-plane aspect ratios above 10. 
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Micron- and submicron- scale patterned ferromagnetic thin films are of interest for appli- 
cations in magnetoelectronicsi""-. Shape anisotropy, through finite aspect ratios of length 
to width, provides a convenient mechanism to adjust the internal field Hd = —NM in 
these structures^, independent of induced anisotropy. In nanometer-dimension patterned 
structures, estimates of the demagnetizing factor from uniform (ellipsoidal) formulae^ are 
considered to represent the anisotropic field well^. 

In uniformly magnetized thin-film structures of 10 fim and greater lateral dimension, ap- 
propriate formulae to estimate demagnetizing fields are not as clear. Two forms for have 
been proposed since the earliest treatments of the topio^"-: "magnetometric" and "fiuxmet- 
ric (ballistic) " , which are the volume average and area average over the midplane normal 
to M, respectively. These two forms diverge for structures with very large aspect ratios of 
thickness b to in-plane dimensions a, c (thin-film limit )^. Nevertheless, the magnetomet- 
ric form is widely applied in the literature^i'^, typically without recognition that alternate 
forms exist. We are not aware of prior attempts to validate either form through comparison 
with the other in fitting experimental data. 

In this paper, we show that the fiuxmetric form is in fact the superior form for the de- 
magnetizing factor N in the thin- film limit. We have evaluated this formula both in the limit 
of saturated magnetization in high field, using variable frequency ferromagnetic resonance 
(FMR), and in low-field hysteresis, using magnetic optic Kerr effect (MOKE) measurement. 
We find agreement of Ho to 20% with a simple analytic limit of the formulae shown in Ref. 
9. We find, on the other hand, that the magnetometric form is a very poor estimate for thin 
films, disagreeing with experimental data by a factor of four for the structures considered in 
this study. 

Patterned Py films (Ni8oFe2o) were deposited on Si substrates using magnetron sputter- 
ing with base pressure better than 3x10"^ Torr. The structure of the films was Si/Si02 
substrate/ Ta(5 nm)/ Ni8oFe2o(40 nm)/ Ta(3 nm). Laser-direct-write photolithography at a 
resolution of 1 /im was used to fabricate rectangular stripe patterns (Fig. la,b), with short 
dimensions ranging from 10 fim to 150 fim. A 1500 fim x 1500 fim square was used as 
the unpatterned comparison. Induced magnetic anisotropy was introduced using an in-situ 
quadrupole electromagnet setup^ in UHV. A rotating field of Hj = 150 Oe was applied 
in phase to the rotating sample holder at 0.25 Hz during sputtering. The Py films were 
deposited together on two identically patterned substrates. The first substrate, denoted as 
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" Py-parallel" (PP), was oriented such that Hi is parallel to the long axes of the elements. 
The second substrate, " Py-orthogonal" (PO), has the orthogonal orientation. The films 
were post-annealed at 250 °C in vacuum of 10~^ Torr for 1 h under a field of 4.0 kOe along 
the deposition field to strengthen the induced anisotropy. 

The hysteresis loops of the stripe arrays on PP and PO were characterized via MOKE 




FIG. 1. (a,b) Patterned Ni8oFe2o stripe arrays with short dimensions ranging from 10 ^m. to 
150 ^m. The distances between adjacent stripes are the same as their widths, (c) Definition of 
Cartesian coordinates and the dimensions. 

with the biasing field applied along the short axes. The size of the laser spot is close to the 
array dimension (1.5 mm). The two substrates were then cut into individual arrays using 
a dicing saw. The in-plane shape anisotropy of the arrays was then characterized by FMR. 
The samples were scanned along their long and short axes. The dispersion curves u){H) of 
each array were acquired with rf frequency varying from 2 Ghz to 26 Ghz. The curves were 
fitted by the Kittel function: 

f = ^ll'li^-nf ■ (ifex + liA){Ile. + Ha + Ms) (l) 

where H^x is the external biasing field, Ha is the anisotropy field, 7/27r = 2.799 Mhz/Oe • 
Qeff/'i, and fioMg is the saturation inductance. The differences of the anisotropic fields 
(AHa) on the two axes are recorded in order to determine the shape anisotropy. 

To estimate Ha and analyze the experimental data, consider a rectangular stripe with 
coordinate system defined in Fig. 1(c), with the long dimension, thickness, and 

short dimension along x, y and z axes, respectively. We present a greatly simplified formula, 
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Nfs, of the fluxmetric demagnetizing factor Nj from the full analytic form (Eq. (2) in Ref. 
9), which could work as a substitute for Nf under thin film limit {b << a, c) and when the 
stripe is uniformly magnetized along z: 

TiN% = — (V4a2 + c2 - c) ^ — (2a - c) (2) 

(XC- etc 

The second approximation is valid when the in-plane aspect ratio is large (a >> c). When 
the magnetization is along x, N^^ is two orders of magnitude smaller than N^^ and does 
not contribute to the shape anisotropy. The full analytic formulae of Nj and from 
Refs. 8 and 9 are compared in Fig. 2 along with A^^^ for the thin film limit and a large 
in-plane aspect ratio (a/c=20). We highlight the significant difference between calculated 
magnetometric and fluxmetric forms Nm and Nf, disagreeing with each other in the analytic 
calculation by a factor greater than three, as plotted in the inset of Fig. 2. We also highlight 
the accuracy of Nf^ to the full fluxmetric form A^^: the filled circles are indistinguishable 
from the line in Fig. 2, disagreeing with each other by a maximum of 0.03% for the range 
considered in the calculation. 

In an array of thin-film elements, additional dipolar fields arise from neighboring ele- 
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FIG. 2. Comparison of Nf, N^ and Nfg in the thin-film limit and at large aspect ratio. Note the 
large divergence between Nf and Nm for wide samples (c > 5 fj,m). Inset: ratio of N^ to Nf. 

ments. The stray dipolar field Hg from the magnetostatic interaction can be calculated by 
summing up all the dipolar fields from the neighboring elements on each element and then 
taking an average of all the local fields. The stray field factor is defined as H^'^ = —M^N^'^, 
with the superscript denoting the magnetization direction. Parts of the calculated ^Nf and 
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ANs are listed in Table 1, with definition of ANj^s = — In the structures ANs is 
about a quarter of ANf but with opposite sign, reducing the total demagnetizing field. 
When the rectangular stripe defined in Fig. 1(c) is under an external field H^^ along z, 

array ANf M^ANf ANs MsANs 

"10/im" 23.4e-4 22.20e -5.8e-4 -5.50e 

"15/im" 15.6e-4 14.80e -4.0e-4 -3.80e 

"30/im" 8.5e-4 S.lOe -2.0e-4 -1.90e 

"150Aim" 1.6e-4 1.50e -0.7e-4 -O.70e 

TABLE I. Lists of calculated ANf and AA'^, for dimensions in Table 1. The equivalent anisotropic 
fields are calculated assuming ^qMs = 0.95 T. 

the total energy can be expressed as: 

E = -fioMsHeccCosO + [±K + ^fioM/{AN* + ANs)]cos^ 9 (3) 

where K is the induced anisotropy energy, 6 is the angle between Mg and H^x, AN* stands 
for either AA^^ or ANf. The sign of the induced anisotropy energy is positive for PP and 
negative for PO. In a classical model, K = fioMsHko/2, where Hko is the induced anisotropy 
field. The new saturation field with shape anisotropy becomes: 

Hk = ±Hko + Ms{AN* + AA^,) (4) 

One should notice that when Hk < the hysteresis loop will show easy-axis behavior. 

In FMR the effective anisotropy field can be calculated by Ha = —OEa/ fJ'odMg where 
Ea stands for the second term in Eq. (3). The magnetization of the elements in an array 
will precess about H^x with small amplitudes. However because the movements of different 
elements are not in phase, the stray dipolar fields are decoupled from the resonance precession 
and we approximate them as a constant field. So in Eq. (3) the term ^fioMg'^ANsCosO'^ 
should be replaced by /loMs^ A']cos 9, where A"* stands for the stray dipolar factor along the 
biasing field direction. The effective anisotropic fields along x and z are calculated to be: 

^A^'' = [-]{±Hko + M.,AN*) - MsN^'^^^ (5) 

At fixed frequency, the difference of the effective anisotropic field between x and z is: 

AHa = H\-Hl = 2[±Hko + M,(AA^* + AA^,/2)] (6) 
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Figure 3 shows the experimental Hk from MOKE, in comparison with the analytic Nf 
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FIG. 3. Comparison of experimental from MOKE with the prediction of Nf (dot) and Nm 
(dash) from Eq. (4) for PP and PO. Sohd hne: data fitting into Eq. (4) by using Nj. Inset: 
Selected MOKE loops of PP arrays. 

(dot) and Nm (dash) from Eq. (4), as a function of 1/c. The experimental Hk of each array 
is obtained by extending the MOKE slopes at if = to saturation. For PO arrays (not 
shown), when the stripe width is larger than 30 yum, cannot be determined. The inset 
shows the hysteresis loops of PO arrays, all of which show linear HA response with He < 
1 Oe. The prediction curves are calculated taking the saturation inductance ^^Mg = 0.948 
T measured in FMR, and induced saturation field Hko = 3.4 Oe measured from the 1500 
/imxlSOO fim square. 

The best fits to the MOKE-derived H^ are shown by the solid lines. Here we treat ^qMs 
and iffco as free parameters but fit the two curves simultaneously. We find that the MOKE 
Hk values are close to the Nf prediction, and far from the Nm prediction. The data closely 
approximate the fitting lines to Nf, validating the simple approximate form in Eq. (4). 
The fitted /loM, is 1.17 T, 23% larger than 0.948 T. The fitted HkQ are 4.5 Oe(PP) and -2.3 
Oe(PO). The former result is close to the 3.4 Oe but the latter is with an opposite sign. This 
is due to the domain wall movement rather than domain rotation and has been discovered 
also in patterned spin valves^^. 

Figure 4 shows the result of FMR measurement. In Fig. 4(a) the f{Hex) relations with 
Hex along X and z are fitted into Eq. (1) in order to extract Qeffi Ha and /ioMg. All the 
arrays can be fitted well by the the same HqMs of 0.948 T and g^ff of 2.115. 
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The differences in Ha are summarized in Fig. 4(b), in comparison with the analytic 
Nf (dot) and Nm (dash) from Eq. (7), as a function of 1/c. The fitting into Eq. (7) is 
also completed with the same ^oMg for the two series as the fitting parameters. Very close 
agreement between the experimental data and fits is seen, with the fitting curves (solid lines) 
very close to the Nf prediction (dot). The fitted H^q are 4.7 Oe for PP and 3.1 Oe for PO, 
both close to the induced saturation field 3.4 Oe. The fitted fioMg is 1.08 T, 14% larger 
than 0.948 T from FMR. 

Summary: we have shown that the fluxmetric demagnetizing factor Nf is superior to 
the magnetometric form A^^ for large- area thin film rectangles (length/thickness > 200). 
Furthermore, we have shown a very simple analytical approximation, A^^^^, to A'^^ in Eq. (2), 
which is excellent over the full range studied. These results will facilitate magnetostatic 
biasing of thin-film structures for applications in magnetoelectronics. 

We acknowledge support from the US Department of Energy grant No. Ide-ee0002892l 
and National Science Foundation No. ECCS-0925829. 
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